Observation of interspecies 6 Li - 133 Cs Feshbach resonances 
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We report on the observation of nineteen interspecies Feshbach resonances in an optically trapped 
ultracold Bose-Fermi mixture of 133 Cs and 6 Li in the two energetically lowest spin states. We assign 
the resonances to s- and p-wave molecular channels by a coupled-channels calculation, resulting in 
an accurate determination of LiCs ground state potentials. Fits of the resonance position based 
on the undressed Asymptotic Bound State model do not provide the same level of accuracy as the 
coupled-channels calculation. Several broad s-wave resonances provide prospects to create fermionic 
LiCs molecules with a large dipole moment via Feshbach association followed by stimulated Raman 
passage. Two of the s-wave resonances overlap with a zero crossing of the Cs scattering length which 
offers prospects for the investigation of polarons in an ultracold Li-Cs mixture. 

PACS numbers: 34.20.Cf, 34.50.-s, 67.60.Fp, 71.38.-k 



The ability to control and manipulate atomic inter- 
actions via magnetically tunable Feshbach resonances 
(FRs) presents a unique toolbox for the field of ultra- 
cold atoms [l|. FRs are used for the creation of deeply 
bound molecules via Feshbach association B 31 followed 
by stimulated Raman adiabatic passage [4|-|6j], which 
gives access to the study of many-body physics, quan- 
tum chemistry and precision measurements 0, H| ■ With 
a permanent electric dipole moment of 5.5 Debye [9l.[l0|. 
the largest among all alkali-metal dimers, a system of 
LiCs molecules in their energetically lowest states [ll[ 
is considered to be an excellent candidate for the in- 
vestigation of dipolar quantum gases 12j]. Another ap- 
plication of the precise tunability close to a FR is the 
study of Efimov trimers The large mass ratio of 

mcs/wLi = 22 results in an advantageous universal scal- 
ing factor of 4.88 instead of 22.7 as found for a system 
of equal masses [l4| , resulting in excellent conditions for 
observing a series of several Efimov resonances, which, so 
far, has not been achieved. Interspecies FR can be also 
used to control the interaction between an impurity and 
a Bosc-Einstcin Condensate (BEC). Such a system can 
directly be mapped to the Frohlich polaron Hamiltonian 
[F5l - [l7| . which describes the interaction of an electron 
gas with the charged lattice atoms in a crystal. As the 
excitations in a Bose-Einstein condensate represent the 
lattice phonons, a FR allows the precise adjustment of 
the modeled phonon-electron coupling strength ev, thus 
allowing one to explore fundamental solid-state systems. 

FRs between different alkalin e sp ecies have been iden- 
tified for a variety of Bose-Bose [18142 1] , Bose-Fermi mix- 
tures [22h28| a s well as for the Fermi-Fermi mixture 
of 6 Li- 40 K 12911 . While the broad intraspecies FRs in 
133 Cs [H Elf and 6 Li [Hll are well-explored, until 
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now, only little is known about the interspecies inter- 
action properties of Li and Cs. Initial theoretical in- 
vestigations of the ground state potential curves by ab 
initio calculations |3, H3| could be improved by spec- 
troscopic data based on Fourier-transform spectroscopy 
!38j and photoassociation spectroscopy [3!| , mainly per- 
formed with 7 Li 133 Cs. The elastic scattering proper- 
ties of 7 Li- 133 Cs were studied by means of thermaliza- 
tion measurements at zero magnetic field [401] ■ however, 
the tunability of the interspecies scattering properties 
via tuning the magnetic field has yet remained unex- 
plored. Here, we report on the observation of nineteen 
interspecies loss features in different spin channels of an 
optically trapped 6 Li - 133 Cs mixture, by scanning a ho- 
mogeneous magnetic field (Feshbach spectroscopy). The 
magnetic field positions and widths of the observed FRs 
are analyzed in a full coupled-channels calculation, al- 
lowing a consistent assignment of the resonances. 

We have realized an all-optical preparation scheme to 
simultaneously trap an ultracold 6 Li - 133 Cs mixture with 
a magnetic field control up to 1300 G by sequentially 
transferring Li and Cs atoms into an optical dipole trap. 
The preparation scheme for fermionic Li follows the ex- 
perimental approach described in 35] . In brief, we trans- 
fer 10 6 Li atoms in a mixture of \f = 1/2, mf = +1/2) 
and 1 1/2, — 1/2) states from a magneto-optical trap 
(MOT) containing 10 s atoms to a crossed optical dipole 
trap (ODT) (1/e 2 waist 60 fim, crossing angle 8.5°, wave- 
length 1070 nm), which creates an effective trap depth of 
Uu = 1.7 mK at a laser power of 150 W per beam. Af- 
ter an evaporation phase at 760 G, where the power of 
the ODT is decreased in a nonlinear ramp to 0.9 W, the 
power of the ODT is again adiabatically increased to 1.8 
W in order to minimize losses of the Li sample during 
the Cs preparation. 

After completion of the Li loading cycle, a Cs MOT 
is loaded slightly displaced to the center of the ODT, in 
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FIG. 1. (Color online) Interspecies 6 Li - 133 Cs Feshbach resonances, (a) Calculated magnetic field dependence of scattering 
lengths for the entrance channels Li 1 1/2, 1/2) ffi Cs|3, 3) (blue line) and Li|l/2, -1/2) ffi Cs|3, 3} (green line), (b) and (c) 
Remaining fraction of Li atoms, normalized to a sample without Cs, in the state 1 1/2, —1/2) and 1 1/2, 1/2), respectively, after 
a hold time of 100 ms at the given magnetic field. The broad loss features at 943 G in b) and 893 G in c) arise from a Feshbach 
resonances of a small amount of Cs in the |3, 2) state, (d) and (e) Narrow interspecies loss features with doublet structures 
in the Li 1 1/2, — 1/2) and Li 1 1/2, 1/2) state, respectively, measured under different conditions (see text) and normalized to 
the offset lines. The data points are an average over five randomized measurements. Center positions and widths of the loss 
features are obtained by fitting a sum of Gaussians (orange lines). 



order to avoid interspecies two-body losses [4l|. During 
a compression phase with increased magnetic field gradi- 
ent and laser detunings as well as reduced laser powers, 
the sample is superimposed with the ODT and further 
cooled and spin polarized by three-dimensional degener- 
ate Raman-sideband-cooling (DRSC) 42j. Reducing the 
ODT power to 0.9 W during the DRSC phase transfers 
the Cs atoms into the dipole trap, while the application 
of an offset field after the DRSC maintains the spin po- 
larization. After 125 ms in the ODT, oscillations of the 
atomic cloud are damped out and the Cs gas is thermal- 
ized. Since the extension of the atomic cloud during the 
DRSC is significantly larger than the trapping volume of 
the ODT, we observe the Cs temperature increase from 
T Cs = 1 (jK after DRSC to T Cs = 8 during the trans- 
fer due to the non optimal spatial mode matching. 

The different polarizabilities of Cs and Li for light at 
1070 nm results in effective trap depths of Uqs/^b = 
40 and U Li /k B = 10 fjK at a power of 0.9 W. We 
observe 5 xlO 4 Cs (2 xlO 5 Li) atoms at temperatures 
Tq s — 8 /iK ( Tu = 2 /iK). The trapping frequen- 
cies were measured via trap oscillations to be wcs/27r 
= (380, 380, 30) Hz and w Li /27r = (900, 900, 65) Hz. 
Prom these measurements we deduce atomic peak densi- 
ties of n Cs = 3 x 10 11 cm" 3 (n u = 1 x 10 12 cm" 3 ). The 
Li sample is spin-polarized by removing one of the two 
rrif components at high magnetic fields with a resonant 
light pulse. The Cs atoms are to 85% spin polarized in 
the |3, +3) state due to the DRSC, with the remaining 
fraction populating the |3, +2) state. 

Feshbach spectroscopy is performed by ramping the 
magnetic field to a desired value and waiting for a vari- 
able hold time between 60 and 2000 ms, which is exper- 



imentally optimized to achieve Li loss signals between 
15% and 55%. For four very weak loss features, the Cs 
density in the ODT was furthermore increased by fac- 
tors of 1.5 — 3 by using longer loading times of the Cs 
MOT. The remaining fraction of Li atoms is monitored 
via high-field absorption imaging providing selectivity on 
the spin state. A FR manifests itself as an enhancement 
of atomic losses due to an increase of the inelastic three- 
body recombination rate [lj . In addition, the increase in 
the elastic scattering cross section on a resonance leads to 
a significant heating of Li by the much hotter Cs atoms, 
resulting in loss of Li atoms because of the lower trap 
depth for this species. Working with comparable Li and 
Cs atom densities, the larger visibility of FRs is thus 
expected by monitoring the Li atom number during the 
measurements. To calibrate our magnetic field we drive a 
microwave (mw) transition of Li at the magnetic fields of 
the maximum of the interspecies losses and determine the 
magnetic field from the measured mw frequencies with 
the Breit-Rabi formula. Fitting Gaussian line profiles to 
the mw spectra yields FWHM of 200 kHz, indicating a 
magnetic field accuracy of 140 mG, probably limited by 
the residual inhomogeneity of the magnetic field over the 
atomic cloud. 

We find a total of nineteen loss maxima for the two 
Li mf entrance channels, all within the range between 
650 G and 950 G. Figure [1] shows four Li loss spectra 
in the field region, where the broadest loss features had 
been obtained, as well as scans over two narrow features, 
showing their doublet structure. The loss features dis- 
appear when the Cs is removed with a resonant light 
pulse during the DRSC or when the measurement is per- 
formed under identical conditions with the other Li mf 
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component. Furthermore, we repeat the measurement 
with only Cs in the trap, to ensure that the observed 
Li losses are not associated with Cs intraspecies FRs, 
which might indirectly influence the number of optically 
trapped Li atoms. Positions and widths of the losses were 
determined by Gaussian fits, and the results are listed in 
Table Q] Trapping both spin components of Li simul- 
taneously with Cs atoms for 500 ms, the full accessible 
field range between G and 1300 G was scanned with 
a step size of 150 mG, yielding no observable additional 
interspecies loss features. 

We analyze the obtained resonances in a coupled chan- 
nels (cc) calculation in a similar way to previous calcula- 
tions on other alkaline systems (see, e.g., [HI Hi]]), using 
the atomic parameters for Li and Cs as compiled by Ari- 
mondo et al. 43]. Precise determination of FRs and 
scattering lengths relies on accurate potential curves for 
the a 3 E + and X x Yj + states of LiCs. Therefore, both 
potentials are constructed in a power series of the inter- 
nuclear separation R, similar to the parametrization of 
NaK ground state curves [Hj]. In a series of iterations 
between cc calculations for the FRs and single channel 
calculations of rovibrational levels within the singlet and 
triplet potential, the description of the potentials is im- 
proved by fitting the underlying coefficients [45[ to match 
the locations of the FRs and simultaneously the 6498 
rovibrational transitions from laser-induced fluorescence 
Fourier-transform spectroscopy of the previous work [38[ . 

The resulting resonances are listed in Table [ I by their 
deviation 8=B**g — B^°°, and the resulting scattering 
length dependence on the magnetic field is depicted in 
Figure [1] B^° is the peak position of the two-body col- 
lision rate at a kinetic energy of 2 /iK. The fit uses the 
experimental uncertainty for the weighting. The fit re- 
sults in a reduced chi-squared of 0.75 for all Feshbach 
resonances, which is very satisfactory, since only 3 res- 
onances are slightly beyond the experimental uncertain- 
ties. We checked that the s-wave resonances will not 
shift by more than —0.02 G, choosing a kinetic energy 
of the entrance channel of 20 nK instead of 2 /iK, as 
it is adequate from the experimental conditions. The 
redefined potential curves reflect scattering lengths for 
6 Li - 133 Cs ( 7 Li - 133 Cs ) of a sing i c t = 30.252(100) a 
(a s in g iet = 45.477(150) ao) in the singlet potential and 
applet = -34.259(200) a (a tri pi c t = 908.6(100) a ) in 
the triplet potential. 

A splitting of four p-wave resonances could be re- 
solved in the experiment, which originates from different 
projection of the rotational angular momentum I onto 
the space fixed axis. This splitting is described by the 
effective spin-spin operator Vdip(-R) for which we used 
the effective form [46]: V dip (R) = |A(i?)(3S"| - S 2 ), 
where Sz is the total electron spin projected onto the 
molecular axis. The function A contains the direct mag- 
netic spin-spin interaction as the first term and sec- 
ond order spin-orbit contribution as the second term 
A = — |a 2 (jly + etso ex P {~bR)), and is given in atomic 
units with a the universal fine structure constant. The 



fit resulted in aso = —8.0 au and b = 0.8 au. 

All FRs originate from the least bound levels below 
the atom pair asymptote 2s 2 Si /2 (/Li = 1/2) + 6s 
2 S , 1 / 2 (/cs = 3) and are significant mixtures of the triplet- 
singlet manifold shown by an expectation value 0.7 of the 
total electron spin. Thus only the projection mp of the 
total angular momentum and the rotational angular mo- 
mentum I are good quantum numbers. The total angular 
momentum /, excluding / at zero magnetic field, and, due 
to the strong Cs hyperfine coupling, the angular momen- 
tum G = S + ics are also good quantum numbers which 
further justifies the assignment we give in Table HI 

To characterize the resonance profiles we give the reso- 
nance parameters for truly elastic scattering channels of 
two body collisions at the chosen kinetic energy: back- 
ground scattering length ab g and the field difference be- 
tween the resonance position and its zero crossing A. 
Both were determined by fitting the conventional func- 
tion of the scattering length as function of the magnetic 
field B: 



using as many terms as resonances are present in this 
channel. The channels Li|l/2, 1/2) © Cs|3,2) and 
Li |l/2, —1/2) © Cs |3, 2} have significant inelastic contri- 
butions, thus such simple representation of the scattering 
length as function of B will not be appropriate. 

Besides two weak and broad features at 937 G and 
988 G in the s-wave entrance channels Li|l/2,±l/2) © 
Cs|3,2) and a very sharp feature at 1019 G for Li 
1/2, 1/2) state, the cc model suggests a series of d-wave 
resonances in the field region below 400 G, which all were 
not detected in the scan over the full field range. No fur- 
ther s- wave resonances for magnetic fields up to 1300 G 
are expected for the studied entrance channels. 

We note that an assignment of the resonances via a 
weighted least-square fit to a simple, undressed Asymp- 
totic Bound State Model (ABM) [Hj, does not pro- 
vide satisfactory agreement with the full set of observed 
resonance positions, reproducing magnetic field values 
of the resonances only within a range of 13 G. The 
ABM yields energies of the / = least bound state of 
eg = 3323 MHz for the singlet and e? = 3945 MHz for the 
triplet state, and for I = 1 ej = 2795 MHz for the singlet 
and e\ — 3455 MHz for the triplet state, and a singlet- 
triplet wavefunction overlap of ( = 0.759. These values 
deviate significantly from the results of the cc calculation 
giving eg = 1566 MHz, e? = 3942 MHz, ej = 1159 MHz, 
e\ = 3372 MHz and ( = 0.866. In order to minimize 
the influence of the coupling between resonant channels 
and continuum of the (broad) s-wave resonances, we fit 
only p-wave resonances, resulting in reduced deviations 
with eg = 2546 MHz, e? = 3978 MHz, ej = 1431 MHz, 
ej = 3626 MHz and C = 0.845. In addition, when using 
the values of the cc calculation in the ABM, we obtain 
the observed resonance spectrum, but globally shifted 
by about -60 G. These large discrepancies between the 
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Entrance channel 


B? c x s p [G] 


AB cxp [G] 


5 [G] 


A [G] 


a bg [a ] 


niF 


f 


G 


1 


Li |l/2, +1/2) 


662.79(1) 


0.10(2) 


-0.04 


- 


- 


7/2 


9/2 


7/2 


1 


© Csj3,+3) 


663.04(1) 


0.17(2) 


-0.02 


- 


- 


9/2 


9/2 


7/2 


1 




713.63(2) 


0.10(3) 


-0.05 


- 


- 


9/2 


7/2 


7/2 


1 




714.07(1) 


0.14(3) 


-0.05 


- 


- 


9/2 


7/2 


7/2 


1 




843.5(4) 


6.4(1) 


0.51 


60.4 


-28.5 


7/2 


9/2 


7/2 







892.87(7)* 


0.4(2) 


-0.11 


4.6 


-28.5 


7/2 


7/2 


7/2 





Li 1/2, -1/2) 


658.21(5) 


0.2(1) 


0.07 


- 


- 


7/2 


9/2 


7/2 


1 


Cs |3,+3) 


708.63(1) 


0.10(2) 


-0.05 


- 


- 


5/2 


7/2 


7/2 


1 




708.88(1) 


0.18(2) 


-0.03 


- 


- 


7/2 


7/2 


7/2 


1 




764.23(1) 


0.07(3) 


-0.06 


- 


- 


7/2 


5/2 


7/2 


1 




764.67(1) 


0.11(3) 


-0.05 


- 


- 


7/2 


5/2 


7/2 


1 




816.24(2) 


0.20(4) 


-0.12 


0.36 


-28.4 


5/2 


9/2 


7/2 







889.2(2) 


5.7(5) 


0.46 


59.9 


-28.4 


5/2 


7/2 


7/2 







943.26(3) 


0.38(7) 


-0.12 


4.3 


-28.4 


5/2 


5/2 


7/2 





Li 1/2, +1/2) 


704.49(3) 


0.35(9) 


0.07 


- 


- 


7/2 


9/2 


7/2 


1 


© Cs |3+2) 


896.6(7) 


10(2) 


0.68 






5/2 


9/2 


7/2 





Li 1/2, -1/2) 


750.06(6) 


0.4(2) 


0.06 






5/2 


7/2 


7/2 


1 


© Cs |3+2) 


853.85(1) 


0.15(3) 


-0.17 






3/2 


9/2 


7/2 







943.5(1.1)** 


15** (3) 


2.21 






3/2 


7/2 


7/2 






TABLE I. List of interspecies 6 Li - 133 Cs Feshbach resonances for combinations of the entrance channels Cs / = 3, mj — +2), 
Cs |3, +3), Li 1 1/2, —1/2) and Li |l/2, +1/2). The experimentally obtained positions and widths are deduced by fitting Gaussian 
functions with the center at the magnetic field B^F an d the FWHM AB cxp to the recorded loss spectra. We deduce statistical 
errors by the 95% confidence interval of the fit. A systematic error for the center position of ±0.1 G originates from the 
calibration plus an additional ±0.1 G from day to day drifts. The cesium density and the hold time were varied in order 
to optimize the lithium loss signal. The resonance positions B^° derived from coupled channels calculations are given as 
deviations 8—B^F — B^° with respect to the observations, showing excellent agreement. The assigned quantum numbers of 
the involved molecular states are explained in the text. For the purely elastic channels of two body collisions at the chosen 
kinetic energy, the background scattering lengths ab g and zero crossings of the scattering lengths with respect to the divergence 
position for 1 nK kinetic energy is shown. (*) We note that the loss feature at 892.87 G overlaps with a Cs 3-body loss minimum 
at 893(1) G [3l]]. This might have a minor effect on the amount of losses due to slightly different 133 Cs densities. (**) The fit 
is not unique, depending on the amount of data points which were considered around the small loss feature at 943.26 G. For 
the cc-calculation the averaged value 942.6(1.0) G from several profile simulations was used. 



cc calculation, precisely reproducing all observed reso- 
nances, and the ABM indicate that an adaption of the 
ABM is required, which is beyond the scope of this paper. 

Three of the observed s-wave resonances provide a 
unique tunability of the 6 Li - 133 Cs mixture. The two 
resonances at 896.6 G and 889.2 G are very close to a 
zero crossing of the Cs scattering length, which originates 
from a very broad s-wave resonance at 787 G [3jJ, |48|, |49j . 
As the fermionic nature of 6 Li suppresses intraspecies col- 
lisions in a spin-polarized sample, an intriguing system 
featuring only tunable interspecies interactions can be 
created. Such a system represents an ideal candidate for 
the study of the Frohlich polaron Hamiltonian [l5l - fl7| 
with a 6 Li impurity in a Cs BEC. The coupling constant 
a scales with the ratio of the impurity-BEC (intra-BEC) 

am ( aBB ) scattering lengths asa~ a iB/ a BB B3- ^elf- 
localized phases of the polarons are predicted when the 
coupling extends a critical value [lg, LL2|; which might 
lead to a collapse of the BEC, if the coupling is fur- 



ther increased [50] . The FR at 843.5 G overlaps with 
the broad intraspecies FR in 6 Li at 832.18 G [36] , thus 
resulting in a system with large scattering length both 
between Cs and the energetically lowest Li state, as well 
as between the two lowest Li states. This scenario might 
be used for the study of self-localization of bosonic impu- 
rities immersed in a two-component superfluid Fermi gas 
[5l| . Additionally, the coincidence between an Efimov 
three-body loss minimum at 893 G (3lj j and two inter- 
species FRs allow for the realization of cooling schemes 
for 6 Li - 133 Cs mixtures where Cs atoms are evaporated 
with minimized losses (3l|, w hile effectively cooling a Li 
sample sympathetically [40]. For reaching a quantum 
degenerate 6 Li - 133 Cs mixture following this approach, 
losses during thermalization have to be avoided. It might 
thus be appropriate to implement species selective opti- 
cal potentials [521 ] using two different wavelengths, which 
would allow for creating similar trap depths for Li and 
Cs. 
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